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<^ • Abstract 

^D , We consider a class of particular solutions to the (2+l)-dimensional nonlinear partial 

I differential equation (PDE) ut + d^^Ux^ — Ux^u = (here n is any integer) reducing it to 

O ■ the ordinary differential equation (ODE). In a simplest case, n = 1, the ODE is solvable in 

.^ . terms of elementary functions. Next choice, n = 2, yields the cnoidal waves for the special 

case of Zakharov-Kuznetsov equation. The proposed method is based on the deformation 

(NJ ' of the characteristic of the equation ut — uux^ = and might also be useful in study the 

higher dimensional PDEs with arbitrary linear part and KdV-type nonlinearity (i.e. the 



00 



% 



nonlinear term is Ux-^u) 



1 Introduction 



Method of characteristics |S] is an effective tool for integrability of first order nonlinear PDEs 

^ ' in arbitrary dimensions. Its modification allowing one to integrate (1-1-1) -dimensional vector 

^ ■ equations (the so called generalized hodograph method) was introduced by Tsarev [1] and 

QO ; developed in refs. [21 [3l H] . A class of matrix equations integrable by the method of characteristics 

^ ■ is proposed in [5]. 

'nI" . In this paper we develop a method of construction of the particular solutions to a new class of 

17^ '. multidimensional PDEs using a deformation of characteristics of already known integrable first- 

T-H ' order equations. We will refer to the following first order (l-l-l)-dimensional partial differential 

> ; equation (PDE) 



0, (1) 

which is a simplest equation implicitly solvable by the method of characteristics. It describes 
the break of the wave profile. It is well known that its general solution is implicitly given by 
the non-differential equation 

x + ut = F{u), (2) 

which is a characteristic of eq.([I]). Here F is an arbitrary function which may be fixed by 
the initial condition. We show that the properly modified equation ([2]) allows one to partially 
integrate the higher-dimensional PDEs, in particular, the (2-(-l)-dimensional PDE 

Ut + dl^Ux^ - Ux^u = 0. (3) 

Note that the term "partially integrable PDEs" has two different meanings. First, the 
PDE is partially integrable if the available solution space is not full (see examples in ref.[6]). 



Second, the PDE is partially integrable if its solution space is described in terms of the lower 
dimensional PDEs [7J. Our nonlinear PDEs are partially integrable in both senses. More 
precisely, we consider a family of solutions to the (2 + l)-dimensional nonlinear PDE (|3]) which 
is described by a nonlinear ordinary differential equation (ODE) with independent variable X2 
where Xi and t appear as parameters. This method may be viewed as a deformation of the 
method of characteristics |B]. 
If n = 1, then eq.dl]) reads 

ut + Uxix2 - wwxi = 0. (4) 

In this case we represent the explicite solution parametrized by one arbitrary function of single 
variable. If n = 2, then eq.( l28p yields a special case of Zakharov-Kuznetsov equation [9l [TOl [11] 

Explicite solutions are represented in terms of elliptic functions (in particular, in terms of the 
cnoidal waves) and are parametrized by two arbitrary functions of single variable. 

General algorithm for (2 + l)-dimensional PDE is represented in the next section, Secj2l 
together with examples of explicite solutions. The asymptotic solutions at large t are derived 
as well. Higher dimensional generalization of the proposed algorithm is discussed in SecJS] 
Conclusions are given in Sec. HI 

2 General form of (2+l)-dimensional PDE reducible to 
ODE 

We proceed with the (2+l)-dimensional nonlinear PDE and prove the following theorem. 
Theorem 1. Let function m be a solution of the following ODE (a deformation of eq.(l2])) 

xi + ut= -V? - dl^u (6) 

(n is an arbitrary integer) in the domain 

ai < Xi < 6i, a2 < X2 < b2, t > (7) 

of the space of variables xi, X2 with some boundary conditions 

B'u = aiiu\^2=a2 + a2iu\x2=b2 + a3jC^X2"U2=a2 + a4i<9^2"U'2=fe = Xh i = 1, . . . , ti - 1, (8) 
where functions Xi ^'^^ constrained by the linear PDEs 

{Xi)t = {Xi)x^t + oiii + a2i, i = l,...,n-l. (9) 

Then the function u is a solution of the nonlinear PDE 

Ut + dl^u^^ - uu^, = 0, (10) 



with boundary conditions ([HI 

Proof. To derive eq. lTTOl) from eq.Q we, first, rewrite eq.([6]) in the form 

xi+ut + w = 0, (11) 



where 

w{u) = d:^u-\u\ (12) 

For the further analysis, we shall mention that eq. (TTO|) admits a representation in the form 

Ut + Wx^{u) = Q. (13) 



Moreover, w satisfies equation 



which may be written as 



wt + d'^w^ - uw^^ = 0, (14) 



(dl-u)iut + w,,) = 0. (15) 

Next, applying the differential operators dt, d^^ and L = d^^d'^ to eq. lTTT]) we obtain: 



E. 



Ex 
El 



u + tut + Wt = 0, (16) 

1 + tu^i +w^^ = 0, (17) 

td:ux,+d:wx,=0. (18) 



Now let us consider the combination 

Et - E,u + El (19) 

which reads 

t{ut + <9" Mxi - uux^) +wt + 9" w^i - uwx^ = 0. (20) 



In virtue of eqs.f lT^ and f lT5|) we have Wt + d'^^^w^ — uw^^ = [d^^ — u){ut + WxJ, so that eq. fl20|) 
can be written in the following form 

it + d:^ - u)^ = 0, (21) 

where il) is defined as follows: 

i) = ut + d'^^Ux^ - uux^. (22) 

Eq. (!2T]) is the ODE for the function ip. However, to get original PDE ( ITOj) . we need the trivial 
solution to eq.( l2T|) : ip = Q. For this purpose we impose the following zero conditions on the 
boundary of domain (|7]) : 

5V = aii^|a.2=a2 + "2*^^1x2=62 + a3i<9*.2^U2=a2 + a4i<9*.2^ 1x2=62 = 0. i = 1, . . . , n - 1. (23) 

Now, using the uniqueness theorem for the linear homogeneous PDE with zero boundary con- 
ditions we conclude that the only solution to the boundary-value problem ( l2T]|23|) is zero, i.e. 
V^ = 0, which is nothing but eq.( 1T0|) . 

We shall note that conditions (!23|l given on the boundary of domain ([7]) constrain the 
boundary conditions for PDE (TTOj) . We derive this constraint assuming that the boundary 
conditions for PDE ( TTOl) are introduced by the operators 5* in eq.( l23ll (i.e. we consider the 



boundary conditions ([H]))- For this purpose we apply operators 5* to E^ given by eq. (IT71) and 
subtract the result from eg. (123!) obtaining 

B\ut - u,,t - 1) = 0. (24) 

From another hand, applying the operator dt — tdx^ to eq. ([8]) we obtain 

B\ut-u.,,t) = {x^)t-t{x^).,. (25) 

Subtracting eq. (!25!) from eq. fl24l) we result in eq.([9]) for the functions Xj(a;i, t), i = 1, . . . , n — 1. 
D 

Below we list several properties of the solution u. 

1. As shown in Theorem 1, the zero boundary conditions (!23|) for the function ip effect the 
boundary conditions for the function m as a solution of the nonlinear PDF (j6]). Namely, 
functions Xi ^^ boundary conditions ([8]) are not arbitrary but satisfy system of linear 
eqs.([nD- This system can be readily integrated yielding 

Xiixi,t) =t{aii + a2i) + Ai{r]), i = l,...,n-l (26) 

ri = t^ + 2xu (27) 

where Ai are arbitrary functions of single argument. 

2. In general, the domain (jTj) of variables xi and X2 might be either bounded or not bounded, 
i.e. ai and 02 might tend to —00 while 61 and 62 might tend to +00. 

3. If the function u is solution to both eqs. llTOl) and ([6]) then it satisfies the linear (1+1)- 
dimensional PDF 

ut = tUx, + 1. (28) 

To derive eq. fl28l) we differentiate eq.([6]) with respect to Xi and subtract the result from 
eq. flTOj) . Note that relation fl28|) on the boundary of domain ([7]) is forced by relations Qj, 
as follows from the proof of Theorem 1, see eqs.( !24]l25l) . 

4. In general, eq.([6]) is a nonlinear ODF. However, it is remarkable that the asymptotic 
M°* of the bounded solution in the domain ([7]) of variables xi and X2 with the bounded 
parameters oi and 02 at large t is described by the linear ODF 



d^^u"' = -u^H, (29) 

1 
ms Xi and 

Xi in definition of rj fl7r|) and write 

7]^' = t\ (30) 



because, in this case, we neglect terms xi and -u^ in eq.([6]). In addition, we may neglect 



Analysis of eq. (l29|) depends on the particular n and on the boundary condition. We 
consider eq. (129!) in Sees J2. II and 12.21 for particular examples of nonlinear PDF (TTOj) with 
n = 1 and n = 2. 



2.1 Simplest case n = 1. Family of explicite solutions of nonlinear 
PDE 



If n = 1, then nonlinear PDE (jTOl) reads 

ut + Wxia'2 - ««xi = 0. (31) 

In turn, eq.(|6]) reads 

Ux:2 = 77^^ — Xi — Ut. (32) 

Boundary condition dH} is represented by a single equation (we take 02 = without loss of 
generality) which we take in the simplest form 

u\^2=o = Xiit,Xi). (33) 

Here Xi has the form given in eq. fl2B]) with an = 1 and 0:21 = 0: 

Xi{xi,t)=t + A,{r]). (34) 



Eq. fl32|) is the first order ODE with separable variables. It may be readily represented in the 
form 



u 



'X2 



2 
and integrated. General solution of eq. (l35ll reads: 



-(m-m+)(m-m„), u± = t±y/rj, (35) 



u{x„x,,t) = ^ + V^ ,_^(^^^,)^..^ - (36) 

Boundary condition ( 133|) in virtue of eq.( l34|) prescribes the dependence of C only on the variable 
77: C{xi,t) = Ciji). This condition reads: 

and may be considered as the definition of Ai{ri) in terms of the arbitrary function C{r]). 
Thus, eq. (15B]) represents a family of solutions of eq.( 1TU]) parametrized by one function of single 
variable: 

u{x,,x,,t) = t + y^Y^CX^- ^''^ 

2.1.1 Asymptotics of bounded solution 

Consider the problem in the domain ([7]) of variables Xi, X2 with the finite bounds ai and 61 at 
large t. Then we shall use eq. (l30|) instead of (1271) for rj. Consequently, the asymptotic solution 
will not depend on xi, which follows from eq. (!38l) for solution u. Let us expand solution (138|1 
in powers of e~^^* and write its first nontrivial term as 

2f 

u = C'^\t)e-^-\ Cit) = -^^y (39) 

5 



where C"''^{t) must be a bounded function of t. The same expression may be obtained directly 
from asymptotic equation (!29|) with n = 1. The boundary condition (!33|) in virtue of eq. (l34|l 
yields 

Ai(t) = -t + C"^(t). (40) 

We see that asymptotic solution (l39|) vanishes at t — i- oo. 



2.2 Case n = 2. Cnoidal waves for a particular case of Zakharov- 
Kuznetsov equation 

If n = 2, then nonlinear PDE (TTOl) reads 

Ut + MziX2X2 - MMxi = 0, (41) 

which is a particular case of the Zakharov-Kuznetsov equation [5] when the transversal part 
of the two-dimensional Laplacian dominates (i.e., if the Debye radios in plasma is small in 
comparison with the th = Cg/wHi, where Cs is the sound velocity and wm is the ion cyclotron 
frequency in plasma p] ) . We choose the boundary conditions ([8]) in the form of two following 
equations 

M|x2=a2 = Xi(^,a;i), u\^^=b2=X2{t,xi), (42) 

where Xi, X2 are given by eqs. fl2^ with an = 0:22 = 1, 0:12 = 0:21 = 0: 

Xi{t,xi)=t + Ai{r^), 1 = 1,2. (43) 

In our case, eq.([6]) reads 

1 



u 



X2X2 



2 
It is equivalent to the following first-order ODE: 



u^ -2tu-2xx). (44) 



Ux2 = ^—/=\/u^ ~ 3tii^ — Qxiu + Ci{xi,t) , (45) 

V 3 ^ 

where Ci is a function to be fixed by the boundary conditions (H^ . Introducing field v related 
with u by the equation 

v = u-t, (46) 

we write eq. (l45|) as 



f^2 = ±^ V4f3 - 12(t2 + 2a;i)w - 4Ci(a;i, t), Ci{xi,t) = 6txi + 2t^ - Ci{xi,t). (47) 
2v3 

The solution of eq. fH7|) is the elliptic Weierstrass function with invariants (72 = 12(t^ + 2xi) and 
^3(Ci(xi,t))=4Ci(xi,t): 

y = u-t = W U^ + C2{x,,t),g2,g3iC^ix,,t))] . (48) 



The functions Cj, i = 1,2 may be found from boundary conditions ( l42l) . Putting successively 
X2 = a2 and X2 = ^2 in eq. lHHj) and using eq. (H3|) we see that Ci must depend only on i]. Thus, 
eq. (l48|) represents solution of eq. (HT|) parametrized by two arbitrary functions of single variable: 

u = t + w(^±^ + C2{v),92,93{CM)y (49) 



Boundary conditions fH2|) in virtue of eqs. fH3|) read 

AM = W (^±^C2{v),92,93{Cr{v))^ , (50) 

AM = W (^±^C2{v),92,93{Ci{v)?j . 

which must be taken as the definitions of the functions Ai{ri) in terms of the arbitrary functions 
C.(r7), z = l,2. 

Now we consider a particular case when all the roots of the polynomial under the square 
root in the rhs of eq.( l47l) . i.e the roots of the polynomial 

v^ -3{t^ + 2xi)v-Ci{xi,t) = 0, (51) 

are real. These roots read 

Pi = d+^, (52) 



P' 2 



P3 



where 



l.,-i - ,73) - i(i - ,v/3), . -a _ ^ (, _ 5) . _,. _ ,„ 

d = 2-'/''(c^ix^,t) + ^Jcf(^i^^~ty^^ . (53) 



Thus, all Pi are parametrized by one function Ci{xi,t) which will be fixed by the boundary 
conditions. These roots are real if |(ip = rj, which is the identity provided Cf{xi,t) — 4r]^ < 0. 
In this case, let 

Vi{xi,t) =mm{pi,p2,P3), V3{xi,t) =max{pi,p2,P3), (54) 

and V2{xi,t) be the remaining root out of the list (^1,^2)^3)5 i-e- "^i < "^^2 ^ "^3- Now we write 
eq.dSD as 

Vx2 = ±^a/(^ - ^i)(^^ - V2){v - V3). (55) 

Let us assume that v is bounded as 

vi < V < V2. (56) 



Then the substitution 



q= I^L^!L^ 0<g<l (57) 

V2 - Vi 



transforms eq.( l55p into the following form 



fe = 4V^V^V/(1-?^)(1-^V), (58) 

k^ = '-^^^. (59) 

V3 -Vi 



Solution to this equation is the Jacoby elliptic function 



^^ /^^-£-^^^_. xj,^^^ (60) 

f 2 - f 1 V 2V3 



Solving eq.f l60|) for u, we obtain 



+ , f \ < f f +^ f +^^ 2 /^ ^ 3^2 Vf 3 (xi , t) - vi {x^ ~ \ 

u = t + vi{r]) + {V2{xi,t) -vi{xi,t))sn ( ± — -= h 62(0:1, t); A; 1. (61) 

Next, boundary conditions (!62|) suggest us to take Ci and C2 as functions of 77, Ci(xi, t) = Ci{ri), 
62(2:1, t) = C2iji)- This means that fj(xi,t) = Viirf), i = 1,2,3. In addition, sn^(x) is an even 
function of argument. Consequently, we may rewrite eq.(l6T]) as (remember that Vi{ri), i = 1,2, 3, 
depend on 61(77) through eqs.( l52]l53l) ) 



^2l^2Vv3{ri)-Vi{ri) 



u = t + vM + {v2{v)-v,{v))sn'{ '^ ^^^ -w^ +c,(^);fc), C72(r7) = ±62(77). (62) 



Boundary conditions fl42|) in virtue of eqs. fl43|) yield 



^i(r?) = vM + (^;2(r7) - vM) sn' (^iV^^S^^M + ^2(77); k) , (63) 



2v^ 



^2(^) = t^i(r/) + (^2(^7) - ^i(r/)) sn (^ — = + 62(77); k 

which must be viewed as the definitions of Ai in terms of the arbitrary functions 6j, i = 1, 2. 
Consequently, eq.(l62l) represents a family of solutions of eq.( l4T|) with two arbitrary functions 
6j, i = 1, 2, of single variable 77. 

2.2.1 Asymptotic of bounded solutions 

We consider the bounded solution on the domain ([7]) of variables Xi and X2 with finite bounds 
tti and 61 at large t. Asymptotic eq.f l29|) reads: 

<l., = -tu""'. (64) 

This equation has the following general solution 

u""' = 6r(xi,t) COS (Vix2 + C^'{xi,t)\ (65) 
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where C°* is an amplitude and Cf^ is a phase. We consider C"^(xi,t) > without loss of 
generality (sign "— " may be embedded in C2*). At large t, we use eq.( l30|) instead of fl27|) for 
rj. Then Ai, i = 1,2, may be considered as functions of t. The boundary functions Xi{t) P3|l 
become bounded if 



Ai{t) = -t + B,{t), i = l,2. 



(66) 



where Bi{t), i = 1, 2, are bounded functions of argument. Boundary conditions (H2l) at t — ?■ 00 
suggest us to take Cf^ as functions of t, so that eq. (!65l) reads 



u 



C1'{t)cos(Vtx2 + C^'{t)). 



(67) 



Thus, asymptotic (1671) does not depend on xi and is parametrized by two arbitrary functions 
of t, C^'^it), i = 1,2. In the asymptotic case, boundary conditions (l42l) in virtue of eqs.(H3ll 
read 

B^{t) = Cl\xi,t) cos (a^Vt + C^' it)), B2it) = C^%xi,t) cos (b2Vi+C^%t)y (68) 

which must be taken as definitions of -Bj(t) in terms of the arbitrary functions Ci{t), i = 1,2. 

Asymptotics of cnoidal waves (1621) . Now we show that asymptotic solution (1671) is nothing 
but the asymptotic of cnoidal wave (l62l) . 

Since rj is given by eq. fl30|) at large t, we neglect Xi-dependence in Cj and take Ci(t) = 



2t^ - tCi{t) with bounded Ci(t) > 0. Then eq.(l53l) yields d"' = t + '-^. Next we write the 

Ci 



roots Pi^ = pi\t^oo, given by eqs. fl52|) . as follows 

^1 



pr = 2t, pr = -t 




P3 



-t+ 




(69) 



Thus, in the asymptotic limit, we have explicite map between f "* and p. 



as 
i ' 



vT = pT 



-t-Vy, vr=pr 




-f, vr=pr = 2t. 



(70) 



unlike the non-asymptotic case, see eqs. fjMJ) . Now we may calculate 



w 2 , 3 1 W g 



'1 



Consequently, for k in eq. fl59|) we have 



^\l 3 



3t+^/f 



— )■ 0, at t — )■ 00. 



(71) 



(72) 



Since sn(a;; k)\k=Q = sin(a;), and sin^x = |(1 — cos2x), we reduce eq. (l62|) to eq. (l67I) where 



cr{t) = -\if, cr{t) = 2C2it). 

In turn, conditions ( l63ll result in eqs.(l68l) provided eqs. (l66 



(73) 



3 Generalization of algorithm to higher dimensions 

The algorithm developed in Secj2]inay be generalized to higher dimension. In this section we 
introduce a large manifold of (M + l)-dimensional nonlinear PDEs with particular solutions 
satisfying either M- or (M — l)-dimensional nonlinear PDEs where the independent variables 
t (and sometimes Xi) appears as parameter. The increase in the dimensionality of the original 
nonlinear PDE is achieved through the increase in the dimensionality of its linear part, while 
the nonlinearity remains the same. 

We shall remark that one could apply this algorithm to larger class of evolution PDEs 
with integro-differential linear part. This would be also important because this structure has, 
in particular, the multidimensional dispersion- less KP. In this respect, we have to remember 
ref.[T2], where another type of particular solutions to this equation were studied and the wave 
breaking for such solutions was explained. However, multidimensional evolution PDEs with 
integro-differential linear part are out of the scope of this paper. 

We formulate the following theorem. 

Theorem 2. Let the differential operator L have the structure L = d^^L, where L{dxi, "^ — 
1,2..., M) is an arbitrary M-dimensional linear differential operator. Let m be a solution of 
the M-dimensional PDE 

x^+ut= -u^ - Lu (74) 

in the M-dimensional domain of the space of the variables Xi, i = 1, . . . , M, 

ai<Xi<bi, i = l,...,M, t>0 (75) 

with the complete set of boundary conditions 

B'u = x^, t = l,...,K, (76) 

where K is some integer and functions Xi? "^ = 1? • • • ? -^; satisfy constraints 

(X.)t - (X,.k = 5*1, i = l,...,K. (77) 

Then the function m is a solution of the (M -|- l)-dimensional PDE 

Ut + Lu — uux^ = 0. (78) 

Proof. To derive eq.f jTHj) from eq.f l7i|) we, first, write eq.f l7i|) as follows: 

xi + ut + w = 0, (79) 

where 

w = Lu- -u^. (80) 

Similar to eq.( !TO|) . eq. (178|) admits a representation in the form 

ut + w^^ = 0, (81) 

which is the M-dimensional PDE. Again, w satisfies the equation 

Wt + Lw — uwxi = 0, (82) 
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which may be written as 

{L - u){ut + w,J = 0. (83) 

Next, we apply operators dt, d^^ and L to eq. flT^ obtaining: 



Et 
Ex 



u + tut + wt = 0, (84) 

1 + tUx^ + Wx^ = 0, (85) 

tLu + Lw = 0. (86) 



El 
Now we consider the following combination 

Et - ExU + El (87) 

which reads 

t{ut + Lu — uux^) + Wi + Lw — uwxi = 0. (88) 



In virtue of eqs. (l82|) and (!83|l we may write Wt + Lw — uwx^ = (L — u){ut + Wxj^)- Then eq. (l88|l 
gets the following form, 

{t + L-u)ip = 0, (89) 

where ip is defined as follows: 

ip = Ut + Lu — uuxi- (90) 

Eq. (l89l) is the M-dimensional PDE for the function ip. Let us impose the zero boundary 
conditions for the function ip 

B'ip = 0, i = l,2,...,K. (91) 

Then, in virtue of the uniqueness theorem for the linear homogeneous PDEs, we conclude that 
^ = 0, which coincides with PDE (ITHj) . 

Boundary conditions (19T!) lead to constraints (!77|) for the boundary functions x-ii i = 
1,. . .,K. These constraints may be derived similarly to the derivation of eqs.([n]). First, we 
apply operators 5* to eq.( 



B^Ex = 0, i = l,2,...,K. (92) 

Subtracting this equation from the eq. &Y\ we obtain 

S«(wt-^«xi-l) = 0, i = l,2,...,K. (93) 

From another hand, applying the operator dt — tdx^ to eq. (!76|) we obtain 

B\ut-Ux,t) = ix^)t-tix^).^, t = l,2,...,K. (94) 

Finally, subtracting eq. (p^ from eq.( p3|l we derive eq.( !77|l for the functions Xi(a;i,t), i = 
..,K. a 
Function u possesses the properties similar to those of the function u in SecJH 
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l,...,K. a 



1. In general, eq.(FMl) is the M-dimensional PDE. However, if the operator L = L{dx2, • • • 5 ^m) 
(i.e. L is the (M— l)-diniensional operator), then eq. flT^ becomes the (M— l)-diniensional 
PDE. Therewith the variable xi appears as a parameter in this equation. 

2. The zero boundary conditions (!9T|) for the function ip effect the boundary conditions 
( ff6ll for the function u through the liner system of (l+l)-dimensional PDEs (ffTj) for the 
functions Xi, ^ = 1, • • • , -f^- 

3. In general, the domain (1751) of variables Xi, i = 1, . . . , M, might be either bounded or not 
bounded, i.e. a, might tend to —00 while bi might tend to cxd. 

4. If the function m is a solution to both eqs. flTHl) and fl74|) then it satisfies the linear (1+1)- 
dimensional PDE fl28|) . To derive this equation we differentiate eq. flT^ with respect to Xi 
and subtract the result from eq. flTS]) . Note that this property on the boundary of domain 
fITS]) is forced by relations fl77|) as follows from the proof of Theorem 2, see eqs. 



5. In general, eq. (l74p is a nonlinear ODE. However, it is remarkable that the asymptotic u"'^ 
of the bounded solution in the domain (1751) of variables Xi, i = 1, . . . , M, with bounded 
parameters oi and 61 at large t is described by the linear PDE 

Lu^' = -u^H, (95) 

because, in this case, we neglect terms Xi and -u^ in eq. (17^ . We also must use expression 



([30]) instead of (EZD for r]. 

Analysis of this equation depends on the particular L and will not be considered in this 
paper. 

4 Conclusions 

We consider a method of construction of the particular solutions to nonlinear PDE f lTOj) . In 
the particular cases, n = 1,2, these solutions are expressed in terms of the either elementary 
functions with one arbitrary function of single variable (n = 1) or elliptic functions with two 
arbitrary functions of single variable (n = 2). The Zakharov-Kuznetsov equation describing the 
ion-sound waves in a low-pressure magnetized plasma is the physically applicable example of 
such systems [9j . We construct the cnoidal waves for the special case of the Zakharov-Kuzntesov 
equation, when the two-dimensional Laplacian becomes the one- dimensional transversal one. 
Asymptotics of these waves in the bounded space-domain are described by a linear PDE yielding 
oscillating behavior. 

Generalization of this method to higher dimensions is briefly discussed. In general, a non- 
linear (M + l)-dimensional PDE from the considered class flTHl) possesses the family of solutions 
satisfying either M- or (M — l)-dimensional differential equation ([71]). Although we consider 
only evolutionary type differential equations, a certain type of the evolutionary equations with 
integro-differential linear part might be also studied in this way. 

Emphasize that we use deformation of characteristic ([2]) in our constructions. However, an 
open problem is whether deformations of characteristics of other multidimensional PDEs might 
be used in a similar way. 

The author thanks Prof. E.A.Kuznetsov for the motivation of this work. This work is 
supported by the Program for Support of Leading Scientific Schools (grant No. 6170.2012.2). 
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